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Abstract: 

In a number of field theoretical models the vacuum angle 9 enters physics in the combina- 
tion 9/N, where N stands generically for the number of colors or flavors, in an apparent 
contradiction with the expected 2n periodicity in 9. We argue that a resolution of this 
puzzle is related to the existence of a number of different 9 dependent sectors in a finite 
volume formulation, which can not be seen in the naive thermodynamic limit V —>■ oo. It 
is shown that, when the limit V — > oo is properly defined, physics is always 27r periodic in 
9 for any integer, and even rational, values of N, with vacuum doubling at certain values 
of 9. We demonstrate this phenomenon in both the multi-flavor Schwinger model with the 
bosonization technique, and four- dimensional gluodynamics with the effective Lagrangian 
method. The proposed mechanism works for an arbitrary gauge group. 



1 Introduction 



Very soon after the discovery of instantons |J in Yang-Mills (YM) theory it has become 
clear || that the latter possesses a hidden parameter 9 whose effects may show up due 
to a non-trivial topological structure of the theory. It is believed that 9 is an angular 
parameter, i.e. physics is periodic in 9 with period 2n. In particular, the values 9 = 
and 9 = 2n correspond to one and the same theory. This is a direct consequence of the 
topological classification of the gauge theories, which is based on the assumption that 
fields are smooth and regular. 

The fact of the existence of this new fundamental constant has immediately posed 
two difficult questions related to the so-called U(l) and strong CP problems. It has been 
argued by 't Hooft that instantons may lead to a resolution of the U(l) problem ||. 
Later, Witten and Veneziano || have found, within the large N c approach, that physics 
should depend on 9 through the combination 9/N c in order for the U(l) problem to 
be solved. On the other hand, a non-zero value of 9 implies |, | a violation of CP 
invariance in strong interactions, which is not observed experimentally. At present there 
is no convincing theoretical argument as to why 9 is so small. Most likely, the strong CP 
problem can not be solved within the strong interaction sector of the Standard Model, 
and will not be discussed here. 

In the present work we address a different, but related, question. As we mentioned 
earlier, a resolution of the U(l) problem suggests that 9 dependence comes in the combi- 
nation 9/N c . To show this, we recall the famous Witten- Veneziano relation || 

/X' = 12 ^F > E vac ~Nlf{9) (1) 

Here f v i is the rf residue and E vac is the energy of the YM vacuum, which is proportional 
to iV| in large N c limit. A resolution of the U(l) problem suggests [fj that f^rnl, ~ 
N c x l/iV c = 0{N° C ). Therefore, the right hand side of Eq.(§ should also be 0{N Q C ). 
This is exactly the case provided we accept that the function / in formula (HD is actually 
a function of the variable 9/N c , rather than of 9 itself. This completes the standard 
argument showing that the 9 dependence should come in the combination 9/N c only. 
There exist many other arguments (based on analyses of 2D CP N ~ X model, SUSY theories, 
etc.) which support the conclusion that the 9 angle always enters physics in a combination 
9/N where N is the number of "colors" or "flavors", depending on the model considered. 

The question we want to raise (and attempt to answer) can be formulated as follows. 
How can one reproduce the 27r periodicity in 9 (which is a strict constraint following from 
the topological classification) for all physical quantities if the same quantities depend 
on 9 through the specific combination 9/N c ? The answer to this question is well known 
in SUSY models Pfl,J5), where it was shown that the 2n periodicity is recovered when a 
discrete number of vacuum states (N for the SU(N) group, iV — 2 for the SO(N) group, 
etc.) is taken into account. The existence of these states in SUSY models is a consequence 
of spontaneous breaking of the discrete chiral symmetry Z 2 n c — > Z 2 (for the SU(N) gauge 
group), which shows up via a formation of the gluino condensate. In this case, the 
different vacua are labeled by the 9 angle as well as a discrete parameter k = 0,1, ...N — 1 
such that the chiral condensate depends on these parameters as (AA) ~ exp ( tg +^ 7rfc ) . 
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Therefore, when 9 varies continuously from to 2ir, N c distinct and disconnected Bloch 
type vacua undergo a cyclic permutation: the first state becomes the second one, and so 
on. All physical quantities are periodic in 9 with periodicity 2tt, as these vacua can be 
just relabeled by the substitution k ^ k — 1 after the shift 9 — > 9 + 2n, keeping physics 
intact. We should note that such a picture is believed to be correct for an arbitrary gauge 
group, irrespective of the existence of the center of the groupQ. 

In non-supersymmetric models such a scenario apparently can not be realized because 
no discrete symmetry which could lead to such degenerate vacua exists in a pure gauge 
theory. The main goal of the present paper is to argue that the pattern of the 9 dependence 
in usual, non-supersymmetric YM theory (gluodynamics) is to some extent reminiscent 
of what happens in SUSY models, although there are important differences between these 
two cases. 

What will be shown is that a discrete number of states, whose presence is crucial for 
the aforementioned mechanism to work, does exist in non-supersymmetric gluodynamics 
when we consider the theory in a large, albeit finite, volume V. These states represent in 
this case local extrema of an effective potential, and have different energies. Under the 
shift 9 — > 9 + 2it, they transform to each other by a cyclic permutation, while some two 
of them cross in energy at certain values of 9 f\. 

Thus, the periodicity in 9 with period 2tt is restored in this finite volume theory. 
However, when the thermodynamic limit V — > oo is performed for a generic value of 9, 
only one state of lowest energy can be seen, as all other states have higher energies and 
therefore drop out in the standard definition 

E vac (9) = - Jim i log Z , 9 fixed (2) 

(This is in drastic contrast with the SUSY case where all N c vacua have the same van- 
ishing energy, and thus all survive the V — > oo limit.) On the other hand, due to the 
superselection rule different states do not communicate to each other (and are absolutely 
stable), and therefore the fact of existence of additional higher energy states could be safely 
neglected, in agreement with Eq.(^), for all physical problems except for one. Namely, 
retaining all these states is necessary for the analysis of periodicity in 9. The values 9 and 
9 + 2tt are physically equivalent for this set of states as a whole. Thus, the fact that the 9 
dependence comes in the combination 9/N c in usual ^ = 00 formulae has nothing to do 
with the problem of periodicity in 9, as those formulae refer to one particular state out 
of this set. As for any fixed 9 the information on all additional states is lost in Eq.(^), 
in what follows this phenomenon will be referred to as a "non-commutativity" of the 
thermodynamic limit V — > 00 and the shift 9 — > 9 + 2tt. Bearing in mind that the states 
of the set do not interact owing to the superselection rule, and to have a correspondence 
with the standard definition (0), in what follows we call the true vacuum state a state 
of lowest energy (for a fixed 9) among this set. When defined in this way, the physical 

1 A related question on a role of the torons |To| , |ll| , , jl3| , |14| |, which are field configurations with 
a fractional topological charge, is not addressed in this paper, see Ref.[jl3| for a list of related problems 
and discussions. 

2 This picture of the 9 dependence is similar to the one advocated by Crewther |l5| , Witten jl6| , and 
Di Vecchia and Veneziano Jl7| for QCD with Nj light flavors. 
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vacuum is periodic in 9 with period 27r, but for different intervals of the values of 9 we 
are talking about a different state from the set as a true vacuum. At certain values of 9 
an exact two-fold degeneracy in this set results in vacuum doubling in the limit V — > oo 
(Dashen phenomenon 



As a warm-up example, we first discuss in Sect. 2 the multi-flavor Schwinger model. 
Using the bosonization approach, we show that the 9 dependence is realized in this model 
in the way just described. An important role of an integer valued Lagrange multiplier 
field, ensuring the quantization of the topological charge, is clarified. 

The rest of the paper is devoted to the problem of 9 dependence in four dimensional 
YM theory. To this end, the knowledge of an infinite series of zero momentum correlation 
functions of the topological density G^G^ is required. We discuss in Sect. 3 methods for 
getting information of this sort by matching short distance and large distance properties 
of the theory (a related discussion can be found in the Appendix). We then construct in 
Sect. 4 an effective Lagrangian (more precisely, effective potential) as the (Legendre trans- 
form of) generating functional for zero momentum correlation functions of the marginal 
operators G^ V G^ U and G^ V G^ V . We shown that an integer valued Lagrange multiplier 
field should be introduced in this effective potential to ensure a single-valuedness and 
boundness from below. The presence of this field imposes global quantization conditions 
on the fields of the effective theory, which reflect the topological charge quantization in 
the original YM action. Sect. 5 deals with the minimization of an "improved" effective 
potential obtained by adding this Lagrangian multiplier field to the theory. We show that 
this procedure yields the above picture of 2n periodicity in 9, in close analogy to what we 
find in the multi-flavor Schwinger model. Conclusions and some discussion are presented 
in final Sect. 6. 



2 6 dependence in TV-flavor Schwinger model 

The main goal of this section is to illustrate most essential technical tools, needed to 
address the problem of 9 dependence in YM theory, on a simple toy model. The multi- 
flavor Schwinger model nicely serves this purpose. It exhibits the above pattern of 9 
dependence and can be analysed in quite a straightforward manner. Though most of the 
results presented in this section are not new, we hope that this discussion may help to 
understand the mechanism ensuring the 27r periodicity in 9 in a more complicated case of 
four dimensional YM theory. 

As is known ||19|| , ||20|| , in the N flavor Schwinger model physics depends on 9 through 
the combination 9/N. At the same time, all physical quantities are periodic functions 
of 9 with period 2ir. We here wish to discuss the way in which these two apparently 
contradictory facts become compatible. In what follows we reproduce some of the results 
of Refs. |1~9[], [p0]| using an approach which emphasizes a key role of an integer- valued 



Lagrange multiplier field and thermodynamic limit procedure in this problem. 

We start with the standard Euclidean action of the bosonized one flavor Schwinger 
model. It takes the form 

r I I i , iF9 

Se = J d M^F 2 + ^<t>? + ^F<p - ^m ? cos(Vi^>) + — ], (3) 
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where F = ^e^F^ and m q is a quark mass. In obtaining this formula we used the 
standard boson-fermion correspondence 

V>7 M V> -> —j=^tivd v ^ , -^cos(\/47T0). (4) 

V 7T 



The key observation is the following. We would like to impose explicitly a global constraint 
ensuring that the topological charge, which is determined by the integral / d 2 3?2^, can 
take only integer values. This constraint can be imposed by introducing an integer-valued 
Lagrange multiplier variable n such that the partition function of the theory is defined 
as the sum over n: 

Z = J2 I DFD(f)e- SE+m27T f d2x £. (5) 

It is clear that we have done nothing wrong by defining the partition function in this 
way, because we introduced only a phase multiplier which always equals 1 for integer 
topological charges. This procedure brings a divergent normalization factor yj 1 which is 
irrelevant anyhow. 

Two remarks are in order. First, the constraint / d 2 x^ — I, I — 0, ±1,±2, ... is 
automatically satisfied due to the identity: 

exp (i2im J d 2 x^j =J2 6 (J ~ l ) ■ ( 6 ) 

As for the second, and most important, remark: the 2ir periodicity in 9 is explicitly seen 
from the general expression (|]). Indeed, a shift in 9 by 27r in Eq.(|3]) can be compensated 
for by a shift in n : n — > n + 1 such that the partition function (|5]) is unchanged. The 
crucial point to make this mechanism work is the definition (^) of the partition function 
Z with the prescription of summing over all n. Such a definition was suggested earlier 



(with quite a different motivation) for the Schwinger model by Smilga [21], and for SUSY 
models (with purposes similar to ours) by Kovner and Shifman 0. 



To study the vacuum structure of the theory, we expand, following Ref . , the fields 
F{x) and 4>(x) in the series over spherical harmonics (a compactification on a manifold 
of volume V is implied) F(x) = J2Fi m Yi m (Q), <p(x) = J2 4>imYi m (£i) and keep only the 
zero Fq mode in what follows. This harmonic F = ^ is fixed by the constraint (§), with 
V being the total volume of the system. Integrating over DF and discarding non-zero 
harmonics of the <fi field, which are irrelevant for the present discussion, we obtain for the 
partition function 

+°° ( 2n 2 l 2 1 
~ o exp | — + il(\^K(f) + 9) + Vm q fi cos (v^^o) j • (7) 



~5~ I 



The 2tt periodicity in 9 is explicitly seen in this representation. However, in order to 
discuss the thermodynamic limit, it is more convenient to use an alternative (dual) rep- 
resentation for the same expression (0): 



-V^ +oo 



2 k 



r.—^ i 1 g 2 V 9 

D<j) Q Y ex p| — — ((j) - + —f=) 2 + Vm q /j, cos 
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where we have used the property of 9^(u, x) function: 



rOO , . , , 2 +00 



6 3 {u,x) = -^= £ €-'= E e- M+2il ™. (9) 



k=— 00 t=— 00 



A generalization of this formula for the case of Nf flavors with equal (and very small, 
m q <C g) masses can be achieved by replacing O — > Nf<p in Eq.(|8|) where, again, we keep 
only the relevant for the vacuum structure part of the partition function: 

+ v5 +00 C 2y _ n ~| 

Z ~ j=k D<Pq S 6XP |-^Wo-^+ ^=) 2 + ym^JV /C os(V4^0 o )| (10) 



Formula (0) was derived earlier [21] in the limit 6 = 0, Nf = 1. 

Now we are ready to discuss the periodic properties of the partition function ( JTTj ) in 
the thermodynamic limit g 2 V — > 00 for the strong coupling regime m g < j. It is clear 
beforehand, without any calculations, that Z is a periodic function of 9 with period 2tt 
for an arbitrary iV/, due to summation over the integers k in Eq. (pTOf) . Now we wish to see 
explicitly how this periodicity works in Eq. (|l~0|) in the limit g 2 V — > 00. For sufficiently 
small 9, only one term k = in the sum ([H]) survives the thermodynamic limit , such 
that 

V^0 O = -— , g 2 V^OO , k = , < # < 7T . (11) 

Nf 



The vacuum energy and topological density condensate for solution (|TTD are 

9 

w f 



E vac (9) ~ -COS(V47T0 O ) ~ "COS 



%F 9 

(-> ~ sin- , 0<tf <7 r. (12) 
When 9 > tt, the next term with /c = 1 corresponds to a lowest energy state, such that 

47T0Q = ^ - 4r , g 2 V -»■ OO , fc = 1 , TT < # < 2tT . (13) 

The vacuum energy and topological density condensate for solution fll3"|) are 

2tt , 



£'m C (6 l ) ~ -COS(V47T0 O ) COs( 

<s> ~ - sin( iv7"iv7 ) '' <,S21 ' < 14 > 

This pattern continues for arbitrary values of 9. For the special case 9 = tt, a cusp 
singularity develops, and we stay with two degenerate vacua which both survive the 
thermodynamic limit (for Nf = 2, this phenomenon was discussed long ago by Coleman 
fjl|). In this case the two vacua are distinguished by the sign of the topological density 
condensate (iF/(2ir)). 
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To summarize, we explicitly see that physics depends on 9 for sufficiently small 9 
through the combination 9/Nf. At the same time, the period of the 9 dependence is 
standard 2ir. This result is in agreement with a very different approach of Ref. ||20||. 

We thus see that the following picture emerges : 
1. The 27r periodicity in 9 holds for an arbitrary value of Nf, even if Nf is a rational, and 
not an integer, number (thought it is not a physical situation for the Schwinger model, 
a similar formula for four-dimensional gluodynamics exhibits such a possibility). The 2n 
periodicity in 9 is always restored due to the summation over integer-valued Lagrange 
multiplier variable n in the partition function (|^). 

2. If we take the thermodynamic limit for sufficiently small 9 in Eq.flTOD from the 
very beginning, we obtain a/47t0o — — w~ once and forever for an arbitrary value of 9. 
Proceeding this way, we would not see other terms with k = 1,2... in Eq. (|10D , which 
are responsible for the restoration of 2ir periodicity in 9, simply because they contribute 
zero to the partition function in the thermodynamic limit. This is exactly what happens 
when one starts with the continuum formulation of the theory from the very beginning. 
In this case all terms with k = 1,2... are automatically discarded. Therefore, the 2ir 
periodicity in 9, as a property of the whole set of solutions, can not be seen in this 
formulation. The thermodynamic limit prescription (|2[) and the shift 9 — > 9 + 2tt are thus 
"non-commutative" . 

3. For each given 9 ^ n ( and m q ^ 0) there is one and only one physical vacuum in 
the thermodynamic limit. By convention (^) we define this vacuum as a state of lowest 
energy among the above set. For 9 = n there are exactly two degenerate states which 
both contribute the partition function in the limit V — > oo. In this sense, at the crossing 
point 9 = 7i the physics is non-analytic in 9. This non-analyticity can not be seen in usual 
V = oo formulae valid for small 9 < n, where a contribution of a lowest energy state only 
is retained. An order parameter which labels two degenerate states at 9 = n is the sign 
of the vacuum expectation value of operator (|^) ~ ±m g sin(-^-). 

4. As we mentioned above, the choice of the physical vacuum as a lowest energy state 
is a matter of convention. Owing to the superselection rule, any state in the discrete 
series (e.g. a next-to-lowest state) could serve as a vacuum as well. However, due to 
the permutational symmetry of the whole set of states under the shift 9^9 + 2tt, any 
redefinition of this sort will result in 27r periodicity in 9 for such a vacuum in the physical 
limit V — > oo. 

In the next sections we will see that the very same picture of the physical 9 dependence 
seems to appear in four dimensional YM theory. 



3 Low energy theorems in gluodynamics 

In this section we start to describe steps which have to be done in YM theory to obtain 



an expression analogous to Eq.(10) for the Schwinger model (the vacuum contribution to 
the partition function or, in other words, the vacuum energy) . This aim requires knowing 
the zero momentum part of the YM partition function, which can not be described within 
perturbation theory. It turns out that this object can be studied using some matching con- 
ditions ensuring consistency of the large distance properties of the theory with its small 
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distance behavior fixed by renormalizability and asymptotic freedom. These matching 
conditions are provided by anomalously broken symmetries through a set of Ward iden- 
tities'] for zero momentum correlation functions of operators describing corresponding 
anomalies. This constitutes what is known as low energy theorems in gluodynamics. Our 
aim in this section is to discuss the low energy theorems in YM theory in order to prepare 
a necessary input for the construction of an effective Lagrangian, which will be carried 
out in Sect. 4. 

In what follows we need two Ward identities for zero momentum correlation functions 
of spin gluon currents in gluodynamics. For the scalar channel case, it was shown long 
ago by Novikov et.al. (NSVZ) that these correlation functions are fixed by renormal- 
izability and conformal anomaly in YM theory. Indeed renormalizability and dimensional 
transmutation ensure that any renormalized zero momentum correlation function of the 
d = D = 4 operator G 2 can only be of the form C (/3(a s ) / (4a s )G 2 ) , where C is a numerical 
constant which depends on the correlation function considered, and the renormalized vac- 
uum expectation value {(3{a s ) / {Aa s )G 2 ) ~ Ay M is the only mass scale in the theory, fixed 
by the conformal anomalyf]. For any given zero momentum correlation function of the 
field G 2 , a value of the particular coefficient C can be found using the dimensional trans- 
mutation formula for renormalized vacuum expectation value of operator O of canonical 
dimension d, written in terms of the bare coupling constant go normalized at the cut-off 
scale Mr. 



(O) = const 



M R exp 



]TT 2 ' 



(15) 



bgl, 

where the one-loop /3-function, (3{a s ) = —ba 2 /(2n) with b = (ll/3)iV c and N c stands 
for the number of colors, has been used. The NSVZ theorem p3[ (with the one-loop 
/^-function) then follows by the differentiation of Eq. ([IB]) , taken for O = —ba s /(8n)G 2 , in 
respect to 1/g 2 . When the full /5-function is retained, it reads |23[ 

lim i J dxe^(0\T{^G 2 (x) ^G 2 (0)}|0) = ~M^G 2 ) , (16) 

Note that the presence of the full /3-function in Eq.flTB]) ensures the renormalization group 
invariance of both sides of this relation. An infinitesimally small momentum transfer 

is introduced in order to select a connected contribution to the correlation function 
(|16|). Eq.flTBD stands for the renormalized correlation function where ultra-violet divergent 
contributions are implied to be regularized and subtracted in both sides of of Eq. fll6|) , see 
the Appendix for a discussion on this point. 

Arbitrary n-point functions of the trace of the energy-momentum tensor 

W) = (^G 2 ) = (^G 2 } + 0(a 2 ) (17) 
4a s o7T 

can be obtained by further differentiating relation fll5|) : 

i n [ dx 1 ...dx n (0\T{a(x 1 )...a(x n )a(0)}\0) = {-A) n (a) , (18) 



3 In SUSY models, a similar use of anomalous Ward identities leads to the well known Veneziano- 
Yankielowicz effective Lagrangian p2[ . 

The fact that terms explicitly containing Ay M do not appear in correlation functions of the operator 
G 2 was checked in Ref.EJ using canonical methods with Pauli-Villars regularization. 
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where, as in Eg. (JTH) , a limiting procedure of the vanishing momentum transfer is 
implied. Note that a regularization scheme in Eqs.([l£j|),(|T8|) is assumed to be the same. 

Let us now address zero momentum correlation functions of the topological density 
operator in gluodynamics. As discussed in detail in the Appendix, the renormalized 
two-point function can be written as 

lim • / dxe^(0\T{^GG(x)^GG(0)}\0)=e(^G 2 ) , (19) 

where £ stands for a generally unknown numerical coefficient (note that its N c dependence 
is expected to be £ ~ iV" 1 , in order to match Witten-Veneziano [§] resolution of the U(l) 
problem). When 9 = 0, Eq.(|19|) is the only possible form compatible with both the 
conformal anomaly and Witten-Veneziano construction. We note that the correlation 
function (|l9|) is defined via the path integral, i.e. with Wick type of the T-product. This 
definition ensures that the nonperturbative gluon condensate in Eqs. (p~6f) , (|T9"|) is the same 
quantity. Perturbative contributions are absent in Eq.(^), see the Appendix. 

In writing Eq.(|T9"D we have assumed that it has the same form, i.e. covariant, for 
any small value of 9 (in the Appendix it is proved only for 6 = 0). Provided this is 
the case, the coefficient £ 2 reiterates, analogously to Eq. fllSP , for all n-point correlation 
functions of GG, as can be seen by the formal differentiation of Eq.(|l9|) in respect to 
9. There are three arguments in favor of correctness of this assumption. First, this 
requirement agrees with the large N c line of reasoning due to Veneziano || where one finds 
that a coefficient standing in the two-point function of the topological density operator 
does reiterate in multi-point correlation functions. Second, this postulated covariance 
of Eq . (|I9|) in respect to 9 goes through a self-consistency check by agreement between 
two different calculations of the 9 dependence of the vacuum energy, when one of them 
is obtained by a straightforward use of Eq. (|l9D (see Eq.(^) below), and another one 
is obtained directly from an effective potential (see Sect. 5). Third, such covariance of 
Eq.(|19p follows automatically with an approach used in f2~| . 

As for the numerical coefficient £ in Eq. (|T9|) , there exist a few proposals to fix its 
value. One of them [2^J is based on the hypothesis of the dominance of self-dual fields in 
the YM vacuum, which suggests £ = 2/6. A different choice, £ = 4/(36), was advocated 
in our work [fj4j] using a one-loop connection between the conformal and axial anomalies 



in the theory with an auxiliary heavy fermion. This line of reasoning was an extension 



of arguments used by Kiihn and Zakharov [2S| to evaluate the proton matrix element 



(p\GG\p). Some further discussion on these matters will be given in Sect. 5. As different 
arguments disagree on what the exact value of £ is, in this paper we prefer to proceed 
with an unspecified parameter £. Fortunately, it turns out that the fact of 2n periodicity 
in 9 can be established without knowing a precise value of £, with the only mild and 
reasonable assumption that the parameter £ is a rational number. In addition, keeping 
an unspecified value of £ in Eq.([B5]) makes it possible to study the 9 dependence for gauge 
groups other than SU(N), at least in principle (see Sect. 5). 

As has been shown in , a combined use of relations flllf ) , flP9] ) enables us to calculate 
the 9 dependence of the vacuum energy and topological density condensate for any number 
of colors N c and small values of the vacuum angle 9 by formal resummations of Taylor 
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expansions in 9 for these objects. The resulting expressions read 

EM = (9\ - ^G 2 \9) = (0| - ^G 2 |0) cos(2£0) = ^cos(2^) . (20) 

OZ7T OZ7T 

for the vacuum energy (here E v stands for the vacuum energy for 9 = ) and 

(6\-?-GG\e) = 2ZE v sm(2Ze). (21) 

87T 

for the topological density condensate (in the particular value £ = 4/(36) was used). 



These formulas seem puzzling as they suggest a "wrong" periodicity in 9 (remember that 
£ ~ A"" 1 ) without any hint at possible singular points 9 ~ ir, which could prevent us from 
making the shift 9 — > 9 + 2n. This might force one to conclude that equations (]20| ) ,(|21| ) 
can not be correct on general grounds, even for small values of 9, and the whole derivation, 
leading to relations (p0|) , (pT|) , was in error. However, as we just saw in the analysis of the 
Schwinger model, a fractional 9 dependence, implied in Eqs. (|20|) , (pl|) , can be in perfect 
agreement with the 2tt periodicity in 9, see Eq.flT2|). What will be argued below is that 
Eqs.(|20|),(21) do not contradict the expected picture of 2ir periodicity of physics in 9 with 
a singular level crossing point at 9 ~ n, for any rational number £ (including, of course, 
both aforementioned choices £ = 2/6 |23| or £ = 4/(36) |24j]). As we have found in the 
study of the Schwinger model, the key to understanding the 9 periodicity problem is the 
analysis of a whole set of disconnected vacuum states. It will be shown in Sect. 5 that an 
accurate transition to the limit V — ► oo while keeping all these states restores the correct 
periodicity and analyticity structure of the 9 dependence, irrespective of a particular value 
£ = any rational number. 



4 Effective Lagrangian for gluodynamics 

The purpose of this section is to construct a low energy effective Lagrangian for gluody- 
namics, which would contain all information provided by the low energy theorems in the 
scalar (|l^) and pseudoscalar fll9f ) channels including all multi-point correlation functions 
of operators G 2 and GG, which can be obtained by differentiating the two-point functions 
(|1|) and (HD, see e.g. Eq.©. 

Before proceeding with the presentation, we would like to pause for a comment on 
the meaning of this effective Lagrangian. As there exist no Goldstone bosons in pure 
YM theory, no Wilsonian effective Lagrangian, which would correspond to integrating 
out heavy modes, can be constructed for gluodynamics. Instead, one speaks in this 
case of an effective Lagrangian as a generating functional for vertex functions of the 
composite fields G 2 and GG. Moreover, only the potential part of this Lagrangian can 
be found as it corresponds to zero momentum n-point functions of G 2 , GG, fixed by the 
low energy theorems. (This effective potential still contains an ambiguity which will play 
an important role in what follows.) The kinetic part is not fixed in this way. Thus, such 
an effective Lagrangian is not very useful for calculating the S-matrix, but is perfectly 
suitable for addressing the vacuum properties^. Specifically, space-time independent fields 

Effective Lagrangians of this kind have been used in supersymmetric theories (see e.g. review papers 
p6|). In particular, the so-called Veneziano-Yankielowicz effective Lagrangian has the meaning just 
described, see 



9 



are amenable to a study within this framework. 

The task of constructing an effective Lagrangian can be considerably simplified by 
going over to linear combinations of original operators^ which enter relations (PS|), (|B3) : 

H -^(- G2 + i ^ Gd ) <*=6^(- G2 -4 G<5 ) ■ <22) 

In terms of these combinations, the low energy theorems for renormalized zero momentum 
Green function, Eqs. ( ]T^| ) and fllPp, take particularly simple forms (for an arbitrary value 
of the vacuum angle 9): 

lim i f dxe iqx (0\T{H(x) H(0)}\0) = -4(H) , 

q^O J 

lim % f dxe iqx (0\T{H(x) H(0)}\0) = -4(H) , (23) 
q->o J 



lim % f dxe iqx (0\T{H(x) H(0)}\0) = . 

9^0 J 



It is easy to check that the decoupling of the fields H and H holds for arbitrary n-point 
functions of H, H. This circumstance makes it particularly convenient to work with fields 

We now wish to construct an effective low energy Lagrangian reproducing at the 
tree level all Ward identities (low energy theorems) for the composite fields H, H, such 
as Eqs.(^) and their n-point generalizationsQ To this end, we consider the generating 
functional of connected Green functions with the space-time independent sources J, J 

expliWU, J)] = y ( DAexp — % - [ dx G 2 + ^ + 2 ™ / dxGG 
„ J 4(r J 32n 2 J 



+iJ J dxH + iJ J dx H 



(24) 



Note the (somewhat unconventional) summation over all integer numbers n in Eq . fl24|) , 
which is analogous to the definition (|5]) for the Schwinger model. This prescription au- 
tomatically ensures the 2ir periodicity in 9 and quantization of the topological charge, 
and is completely equivalent to the way the vacuum angle 9 has initially appeared in YM 
theory [0 . The above form of introducing the 9 angle in the path integral will help us to 
understand how the 9 parameter should be installed in the effective Lagrangian formalism. 

We next define the effective zero momentum fields (here and in what follows J dx = V 
is a total 4-volume) 

dW r - dW 
dxh = ~dJ ' J dxh = -qj ' ( 25 ) 

satisfying the equations 

J dxh = (J dxH) , j dxh= (J dx H) . (26) 

6 In this section we change the normalization of the gluon field in comparison to that used in Sect. 3 
by the rescaling — * (l/g)A fl , and use the one-loop /3-function. 

7 For the case of one real "dilaton" fields a = — 6a s /(87r) G 2 , a similar problem of constructing an 
effective Lagrangian was solved long ago by Schechter |^7| , and Migdal and Shifman (see also pj| ) . 
Our derivation below is akin to the one suggested by Cornwall and Soni in p9|. 
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The effective action T(h, h) is now introduced as the Legendre transform of the generating 
functional W(J, J): 

T(h, h) = -W(J, J) + J dxJh + J dx Jh (27) 

which implies 

dT > > (28) 



d J dxh d J dxh 

^From the definition Q2"^) ) and the low energy theorems (|23|) (and their extensions for 
arbitrary n-point functions, see Eq. fll8|) ) we obtain 

Qn+l 



W 



j=j=o 



fdxdx 1 ... dx n (T {H{ Xl ) . . . H{x n )H{0)}) = (-4) n J dx(H) 



d jn+l 

^-^W\ J= j =0 = i n j dxdx 1 ...dx n (T{H(x 1 )...H(x n )H(0)}} = (-A) n j dx(H) 
W\j=j=o = (29) 



dJ k dJ l 



(as before, the connected parts of the Green functions are implied in Eq. (P9|) ). These 
equations are solved by the function 

W(J, J) = -\J d ^ (H)e- 4J - i y 'dx (H)e- 4J . (30) 
Using Eq.(^), we can express the sources J, J in terms of the fields h, h : 

j -Hif)- J -Hif)- (31) 



Inserting these expressions back to Eq . (p0|) , we obtain W as a function of the fields h, h. 
Now the definition (|2"7|) turns into the differential equation for the effective potential 
U(h,h) = -(l/V)T(h,h) : 

This equation is a complex extension of a real differential equation for the "dilaton" effec- 
tive potential of Refs. $I7\ , . Here comes the aforementioned ambiguity of the effective 
potential. Let us compare Eq.(p2|) with the equation for the real "dilaton" field of [27 , 



dU 1 

"W-to—f. (33) 

Eq.(^3|) has the only solution U(a) = (l/4)cr(log a + const). It is the appearance of the 
multi-branched logarithmic function of a complex argument in Eq.(|3"2"D that gives rise to 
the ambiguity which was absent in the real equation d33|). Let us analyse the way it 
appears when Eq. (|32"D is solved. One obvious solution of Eq.(|32"D is 

U^h, h) = h\og^ + h\og^ + D(h-h), (34) 
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where C, C, D are arbitrary complex constants which may depend on (H), (H), 9. How- 
ever, fl3~4|) is not a single- valued function, and is not bounded from below. After the phase 
rotation h — > hexp(27rir) with an arbitrary integer r, the potential ([53]) transforms as 

— — 77TT — 

U 1 (h,h)^U 1 (h,h) + —(h-h) (35) 

which is physically unacceptable. 

A way out in this situation is to sum over all integers r in the partition function, as 
was suggested by Kovner and Shifman [g] in a similar problem arising with Veneziano- 



Yankielowicz effective Lagrangian |22] for SUSY gluodynamics. Yet, in our case this is 
not the end of the story. Indeed we find that there exists another possible solution of 
Eq.© : 

U 2 {h,h) = ±h]og^£j +^!°g(^) +D(h-h), (36) 

where a is an arbitrary real number. From now on we concentrate on the case when a 
is a positive rational number, a = p/q, where the integers p and q are relatively prime. 
Using the formula 

log z p / q = - Log z + 2ni(n + k~) , n = 0, ±1, . . . ; k = 0, 1, . . . , q - 1 (37) 
q q 

( here Log stands for the principal branch of the logarithm), we see that the second form 
(j36|) makes no difference in comparison with fl34|) when only the principal value of the 
logarithm is considered. However, the theories, described by the effective potentials (|34| ) 
and ([36]), are different quantum mechanically as they imply different rules of a global 
quantization for the fields h, h. This quantization arises when the single-valuedness of 
the partition function is ensured by a summations over the integers. As will be shown 
in the next section, it is the second choice ([36]) for the effective potential that can be 
made consistent with both the 9/N dependence and 2tt periodicity in 9 when a proper 
treatment to global quantization constraints and the thermodynamic limit is given. 
We therefore consider the function 

U 3 (h,h) = jhLog^ + ^hLog^ + j(k + n^j (h-h) , (38) 

which satisfies Eq . (|32"D (i.e. the Ward identities (|23|) and their n-point generalizations) 
for any values of the integers n, k from the range n = 0, ±1 . . . ; k = 0, 1, . . . , q — 1. The 
last term in (]38| ) is a particular form of the last term ~ D in Eq. (|36|) . It can be seen that 
arbitrary values of the coefficient D would be uncompatible with the quantization rules 
imposed by the summation over the integers n, k in the partition function. 

Finally, we have to figure out how the 9 angle should be installed in the effective 
potential (0). An answer to this question can be deduced by comparing with Eq. (|24|) . 
In the YM partition function, the 9 angle enters in the combination 9 + 27m, while the 
summation over the integers n is necessary because of a multi-valuedness of the YM 
action in respect to large gauge transformations. Analogously, the last term in Eq.(^) is 
the only one that can accommodate the 9 parameter in the same combination 9 + 27m. 
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The summation over the integers n is enforced this time by the multi-valuedness of the 
logarithm in Eq. ([38|) . As the presence of the 9 angle is implicit in the constants C, C, we 
now make it explicit in the above way and finally obtain the (Minkowsky space) improved 
effective potential F(h, h) by the summation over n, k in the partition function : 

e -iVF{h,h) _ 

+ 

where the constants C, C' are independent of 9 and can be taken real, C — C' = 2eE, 
where E is some positive constant. The improved effective potential F(h, h) is consistent 
with all constraints imposed by the low energy theorems and, by construction, is a single 
valued function possessing the 2tt periodicity in 9, which was present in the initial YM 
partition function. As is seen from Eq.QBTf), the structure of the effective potential F 
is such that it contains both the "dynamical" and "topological" parts (the first and the 
second terms in the exponent, respectively). We would like to note that Eq.(|3"9"D is a direct 
analog of a similar construction for SUSY models |2^, [J. Namely, the "dynamical" part 
of the effective potential is rather similar to Veneziano-Yankielowicz (VY) [g2| potential 
~ M 2 / 3 logM, while the "topological" part is analogous to an improvement of the VY 
effective potential, suggested by Kovner and Shifman ||. We stress that the improved 
effective potential ( p9|) contains more information in comparison to that present in the 
Ward identities (|3|) just due to the appearance of this "topological" part in Eq . (|39"D . 
Without this term Eq.(|3Tf) would merely be a kinematical reformulation of the content 
of the Ward identities (|23|). As will be shown in the next section, this improvement of 
the effective potential turns out crucial for unravelling the correct periodicity in 9 in YM 
theory. 



V? V- / iV (h T h 

n=-oofc=0 I \ 



c> 



h Log 



h 



T /, q 9 + 2ixn\ h-h 
mV \k + - 



(39) 



5 Minimization of effective potential 

In this section a ground state of the dual low energy theory will be determined by a 
minimization of the improved effective potential (IEP) F(h,h) given by Eq. (|3"9"D . Our 
purpose is to find the 9 dependence of the vacuum energy which is defined as a minimum 
of IEP F(h, h). In this calculation the total space-time 4- volume will be kept finite, while 
a transition to the thermodynamic limit V — > oo will be performed at the very end. 
We start with introducing the "physical" real fields p, i] defined by the relations 

h = 2Ee p+ir > , h = 2Ee p ~ iv . (40) 

(This definition implies F{j] + 2it) = F(i]). As will be seen, this condition of single- 
valuedness of the rj field is satisfied with the substitution (|40|).) In these variables, the 
"dynamical" part of Eq.(^) can be written as follows: 

IhLog^ + hLog^J = -iVEe p [(p-l)cos7]-r]smr]] . (41) 
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The summation over the integers n in Eq.fl3D|) enforces the quantization rule due to the 
Poisson formula 



X] M ~ VEe p sin rj — m 

m \P j 



(42) 



which reflects the quantization of the topological charge in the original theory. Therefore, 
when the constraint (42) is imposed, Eq.(|4lD can be written as 



iV (x. r h T r h \ 

hLoq h h Loq - 

4 V y 2eE y 2eEJ 

Using (@2|),(||), we put Eq.© in the form 



-iVF 



-oo q—l 



y~] ^2 5(VE - e p sin rj — m) exp 



m=-oo fc=0 

where we denoted 



-iVE e p (p — 1) cos 77 + im - rt 

q 



-iVE e p (p — 1) cos rj + imlOk 



(43) 



P 



9 k = 9 + 2tt - k . 

q 



(44) 
(45) 



To resolve the constraint imposed by the presence of 5-function in Eqs.(^2|),(44), we 
introduce the new field M by the formula 



5{VE - e p sin rj - m) oc f DM exp ^iMVE e p sin r]-iM-mj 
Going over to Euclidean spaceQ by the substitution iV — * V, we obtain from Eqs. 

^ ( +oo q-l 

F(p,r},M) = -— log I J2 ex Pl~ VEeP i(P~ l)cosr/ - Msinr/} 

I m=—oo k=0 



(46) 



+ im 1 9 + 27ik- + - r] - - M) - e 

\ q q q J ve 



(47) 



Here we introduced the last term to regularize the infinite sum over the integers m. The 
limit e — > will be carried out at the end, but before taking the thermodynamic limit 

V — » oo. Note that Eq. (|47|) satisfies the condition F(j] + 2n) = F{rj) which should hold 
as long as rj is an angle variable. We also note that the periodicity in 9 with period 27r is 
explicit in Eq . (f4~7D . 

Proceeding as was done for the Schwinger model, to discuss the thermodynamic limit 

V — > oo we use the identity @ and transform Eq.(f4~7D into its dual form 

I ( +oo q-l 

F(p,rj,M) = -— log < Yl ex P [-VEe p {(p - 1) cos rj - Msinr/} 

^ ln=-oofc=0 



VE 

IF 



9 + 2nk - + -ri — - M — 2nn 

q q q > 



(4£ 



3 This is not really necessary. All formulas below can be worked out in Minkowsky space as well. 
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where we have omitted an irrelevant infinite factor ~ e~ x I 2 in front of the sum. Eq.( 
is the final form of the improved effective potential F, which represents the YM analog of 
Eq. ( |TUD for the Schwinger model. To discuss the vacuum properties, the function F should 
be minimized in respect to the three variables p, 77 and M. In spite of the frightening form 
of this function, its extrema can be readily found using the following simple trick. As 
at the extremum points all partial derivatives of the function F vanish, we first consider 
their linear combination in which the sum over n, k cancels out. We thus arrive at the 
equations 

OF 

— — = Ee p (pcosr] — M sin 77) = 
op 

dF dF 

— + —— = -Ee p (psmr} + Mcosr}) = 0, (49) 
or/ oM 

which is equivalent to p 2 + M 2 = 0. Therefore, these equations have the only solution 

(p) = , (M) = , (50) 

while the minimum value of the angular field 77 is left arbitrary by them. The latter can 
now be found from either of the constraints dF/drj = or OF /DM = 0, which become 
identical for (p) = (M) = 0. The resulting equation readsf] 

Y y (e + 27rk- -2nn+- rj + 2e- sin ri) 

r^oc^oV 1 Q P J 

x exp jvEcosr? - ^ (e + 2nk- + -ri-2irnj | =0 ' ( 51 ) 

in which we have to take the limit e — » at a fixed 4-volume V. 

One can see that non-trivial solutions of Eq. ([H]) at e — > are given by 

(V)i = --0 + — / + 2vrr , I = 0, 1, . . . ,p - 1 ; r = 0,±l,... (52) 
p p 

Eq.(|52|) shows that there are p physically distinct solutions, while the series over the 
integers r in Eq.(|52"D simply reflects the angular character of the r] variable, and is thus 
irrelevant. By the substitution of Eq. (|52D back to Eq.(fP|) we obtain the energy spectrum 
for the finite volume theory: 

Ei = F(p = M = 0,rj= (rj)i) = -E cos(r])i = -E cos (--9 + — 1) . (53) 



Thus, we have found that the improved effective potential fl48|) has not one, but rather p 
physically different local extrema, when we look at the theory in the finite volume. The 

9 Note a remarkable similarity between Eq.(]5l|) and the equation /i| sm 4>i = (a/N)(6 — 52 "A?) (where 
/ii and 4>i are the masses and phases of goldstone fields, respectively), obtained by Witten JL6| as a 
minimization condition for the effective chiral Lagrangian for QCD. The limit s — > in Eq. 
analogous to the chiral limit — > in this equation. 



1 is 
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fact that the states (|52p have different energies for generic values of 9 is very important. 
This is where we find an essential difference of non-super symmetric gluodynamics from 
its supersymmetric extension. In the latter case, there are N c degenerate states which all 
survive the infinite volume limit, and correspond to the physical Z2N C symmetry of SUSY 
YM theory [^, |9|, . The absence of degeneracy between the states ( |52"D is therefore very 



natural, as there are no discrete symmetries for non-supersymmetric gluodynamics in the 
thermodynamic limit V — > oo where we should stay along with just one true vacuum. 
Yet, as we will see in a moment, retaining the whole set of local extrema (|52|) is important 
to recover the correct periodicity in 9 in the limit V —>■ oo. 

The remarkable fact about the extrema fl5"2|) is that they are related to each other by 
a cyclic permutation under the shift 9 — > 9 + 2n. The physics is perfectly periodic in 9 
with period 2tt, as the minima (t])i, interchanging under the shift 9 — > 9 + 2tt, can be 
just re-labeled without altering anything. One of the minima always has a lowest energy. 
For example, if < 9 < n/q, it is the I = solution in Eq.flS^D- At the same time, we 
observe level crossing with a two-fold degeneracy at certain values of 9. One series of the 
level crossing points is given by 9 = n (mod2ir) ) irrespective of the values of the integers 
p, q. For example, at the first point 9 = it in this series, the I = and I = q solutions 
have the same energy —Ecos(irq/p). This is the same series of level crossing points as 
was found for the Schwinger model. The difference from the Schwinger model is that now 
these values of 9 do not correspond to level crossing of lowest energy states among the 
set (0). Instead, this happens for another series of level crossing points in 9, which is 
different from the previous one as long as q ^ 1. As can be seen from Eq. fl53|) , it is the 
points 9k = (2k + l)7f/q , k = 0, 1, ... ,p — 1 where the lowest energy state is changing 
from the kth to the (k + l)th branch in the set ([52]). 

Let us now see what happens when the thermodynamic limit V — > oo is taken. The key 
observation is that a lowest energy state, which is the only one that should be retained 
in the limit V — > oo according to our convention (fj), corresponds to different minima 
from the set ([52"D, depending on an interval of variation of the vacuum 9 angle. Thus, to 
perform the thermodynamic limit, we should first fix an interval of 9 (say, < 9 < ir/q), 
and only then select the state of lowest energy among the set (|5"2"|). This solution will be 
the one corresponding to the single vacuum state in the limit V — > oo, for all values of 9 
from this interval. This procedure can be described by the formula 



1 (j2 ie xv\VEcos(-l9 + fl) 
hm — lor ' 



v^oo V I 1 + 5, 



0,l,...,p-l (54) 



(here ^-vr/q.o is the Kronecker symbol, equal to 1 if 9 = (2k + l)7i/q or otherwise). 
The multiplier 1/(1 + (^-Tr/g.o) accounts for the two-fold vacuum degeneracy at the points 
9 = (2k + l)ir/q , k = 0, 1, . . . ,p — 1. We note that Eq.([54"D is perfectly periodic in 9 with 
period 2n. 

Eq.(|54]) shows that in the limit V — > oo cusp singularities occur at the values 9 = 
(2k + l)7r/g, where the lowest energy vacuum state switches from one analytic branch 
to another one, much as it occurs in the Schwinger model. One should note that there 
is no physical jump at 9 = n/q. It is rather re-labelling of a lowest energy state. The 
first derivative of the vacuum energy, which is proportional to the topological density 
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condensate, is two- valued at these points. This means that whenever 9 = (2k + l)n/q, we 
stay with two degenerate vacua in the thermodynamic limit (Dashen phenomenon fl8|| . 
see below). This picture of the singularity structure in 9 resembles the one found for the 
lattice Z p model in 4D |30 . 



If, on the other hand, the thermodynamic limit is performed for a fixed value of 9, 
any information on other states is completely lost in Eq.(|54|). Correspondingly, the 2n 
periodicity in 9 is also lost in infinite volume formulae. We have no chance to know about 
additional states when we work in the infinite volume limit from the very beginning. As a 
result, usual V = oo formulae become blind to the very existence of a whole set of different 
vacua, which is just responsible for restoration of the 27r periodicity in 9. Instead, formulae 
corresponding to the formal limit V = oo look as suggesting a "wrong" (different from 
2tt) periodicity in 9 , see e.g. Eqs.(^,^T|). Now we know that this procedure of the shift 
9 — > 9 + 2tt in the U = oo formulae is simply misleading as it is equivalent to going along 
a single analytic solution of the minimization equation (p>T|), which does not corresponds 
to a lowest energy state for a shifted value of 9. Comparing Eqs.(pO[) and (|53|) , we see 
that the former may well describe the 9 dependence in the physical limit V — > 00 for 
small values 9 < n/q. To this end, we should set the ratio q/p, which so far was arbitrary, 
to the value q/p = 2£ (and take E = —E v ). At the same time, analyticity in 9 of each 
separate branch (|5^) shows that the procedure of a formal re-summation of the infinite 
Taylor series for small 9 < rr/q, which has led to Eq . (|20|) [p4j] , is legitimate. We therefore 
conclude that Eq. fl2"0|) , which should be understood as standing for 9 < ir/q, is not in 
conflict with general principles of periodicity and analyticity in 9 for any rational value 
of the parameter £, including both particular choices £ = 2/6 ]23[] or £ = 4/(36) [p4fl . 

Although the problem of fixing the correct value of the parameter £ is beyond the scope 
of this paper, we can not refrain from pausing for a few comments on these matters^. The 
parameter £ is related to a number of different sectors of the theory, which are disconnected 
due to the superselection rule. One could think that this number of sectors is proportional 
to 6, as the formulae £ = 4/(36) |24[] or £ = 2/6 suggest. However, the analysis of 
SUSY theories shows that it might not literally be the case. If in SQCD we change a 
number of flavors Nf keeping a number of colors N c unchanged, the number of sector 
remains the same and strictly equals N c (and not just proportional to N c ), though the j3- 
function 6 ~ 3N C — Nf changes. As a result, the angle 9 enters physics in the combination 
9/N c for arbitrary Nf. Another argument comes from the analysis of softly broken SUSY 
gluodynamics where the gluino is given a small mass m. As discussed by Shifman f26fl , 
the situation is under control as long as m is small, and the number of different sectors 
remains the same N c (though the degeneracy between them is lifted). On the other hand, 
the case of pure YM theory corresponds to the limit m — > 00 (which means physically 
m ^> A). If the number of sectors N c for small m does not discontinuously changes when 
m becomes large, m ~ A, we end up with different sectors for usual YM theory. 

To summarize, different lines of reasoning lead to different values of £ (though they 
all imply £ ~ A"" 1 ), and correspondingly to different values for a number of sectors in 
pure YM theory. However, irrespective of this particular number (we only assume it be 
a rational), we know that there is only one true vacuum in the thermodynamic limit 



10 The arguments discussed below are due to A. Vainshtein ]3l| ] , to whom we are indebted for sharing 
with us his insight. 
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V — > oo and that the period of the 9 dependence is always 2n. 

After this digression we wish to discuss relations for the topological density conden- 
sate. Namely, we would like to see whether Eq.flSTD, obtained by a direct evaluation of 
correlation functions in YM theory, is consistent with results of this section. It is easy to 
see that it is indeed the case. Differentiating Eq.flSHD in respect to 9, we obtain 




-?§ = q -Esm(-l9 + -l) . (^) 
39 p \ p • 

As only the I = term survives the V — > oo limit for < 9 < n/q, we obtain, in agreement 
with Eq.(f2l|) 

1 .(GG) = -^Esm(^-9) . (56) 



327T 2 ' p \p 

Similarly to what occurs in the Schwinger model, for the special case 9 = ir/q we stay in 
the limit V — > oo with two degenerate vacua which are distinguished by the sign of the 
topological density condensate: 

|1> = \Vi=o = -~) , -^(GG) 1 = -^Esinf-) 
p S2tc z p \p J 

|2) = h=i = ~) , ^~ 2 {GG) 2 = + q -E sin (-) (57) 
p 32n 2 p \p J 

As a CP transformation reverses the sign of 9, it exchanges the vacua |1) and |2): 
CP\1) = |2) , CP\2) = |1) . Therefore, the CP symmetry is broken at 9 = n/q. A 
similar phenomenon of vacuum doubling occurs for any point of the form 9k = kn/q, 
k = 1, 2, 3, . . .. For example, the / = and I = q states are analogously related by a CP 
transformation at 9 = ir. The reason we concentrate on the level crossing point 9 = n/q 
is that at this value of 9 the true vacuum (lowest energy state) switches from the I = 
to the I = 1 branch, while at 9 = ir some excited states cross in energy. Therefore, in 
the sense of the lowest energy state among the set (|52"D, the value 9 = n corresponds to a 
regular, CP conserving point. 

Finally, we would like to briefly describe the case of YM theory with an arbitrary 
(orthogonal, exceptional etc.) gauge group G instead of the unitary SU(N) that was 
discussed so far. For such a gauge group the second Casimir constant C^G) and the 
/3-function would be different. Therefore, the only difference from the previous analysis 
in this case would be different values of the integers p and q, while the pattern of 9 
dependence would remain the same. Thus, the mechanism suggested in this paper seems 
to be valid for any gauge group. 



6 Conclusions 

The most important results of the present analysis are the following: 
1. We have demonstrated that physics is periodic in 9 with period 2ir for an arbitrary 
gauge group. This behavior follows from our definition of the partition function for both 
the original and effective Lagrangian formulations of the theory, where the summation over 
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all branches of a multivalued action is imposed. In the effective Lagrangian framework, 
this prescription is necessary for a single-valuedness and boundness from below of an 
effective potential. 

2. The periodicity in 9 with period 2ir is perfectly compatible with the 9/N dependence 
found in a number of models at small 9. The correct periodicity in 9 is recovered when 
a whole set of different branches is taken into account. The standard definition of the 
thermodynamic limit selects only a lowest energy state among this set. As a result, the 
thermodynamic limit and the shift 9 — > 9 + 2n do not commute (in the sense explained 
in the Introduction). 

3. For generic values of 9, there is one and only one vacuum state in the thermodynamic 
limit. For 9 = n/q there are exactly two degenerate states, which are distinguished by 
the sign of the topological density condensate. 

We would like to end up with some speculations. We emphasize again that 9 = n/q 
is a very special point because of the vacuum degeneracy which does not follow from any 
obvious symmetry of the original Lagrangian. This degeneracy may imply the existence 
of domain walls in the theory at 9 = n/q, which are static field configurations depending 
only on one spatial coordinate. An effective potential describing domain walls could be 
obtained from Eq.(48) by freezing the p and M fields. Such a potential is a complicated 
function of the rj field which, however, reduces to the standard Sine-Gordon form near the 
points (p2|), see Eq.(|53]). Provided a kinetic term (we recall that kinetic terms are not fixed 
by the Ward identities) is added, the theory could sustain domain wall configurations. 
The existence of these solutions in gauge theories could have interesting consequences for 
cosmology and particle physics. 

We would also like to speculate that the above "p" non-equivalent states, could be 
really observed in some nonequilibrium high energy processes with a finite geometry, where 
the superselection rule can not be applied (a similar comment was made by Shifman [EB] ) . 
Such a situation could be realized e.g. in nuclear-nuclear collisions, where the appearance 
of droplets of a "false vacuum" would be similar to the production of droplets of disoriented 



chiral condensate, see e.g. [32 for a review. Instead of an arbitrary direction of the chiral 



condensate for the latter, in the former case we would deal with "p" different values of 
the topological density condensate (GG). One expects that this phenomenon, if exists, 
should be related to the physics of the rf meson and CP violation. Yet, it is not known 
at the moment how to formulate this problem in an appropriate way. 

The inclusion of the light quarks into the effective Lagrangian framework and the 
resulting picture of the 9 dependence in QCD will be discussed in [53] . 
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Appendix 



The purpose of this appendix is to discuss in somewhat more detail the derivation of 
the low energy theorems (|l^) and (|T9|) and, in particular, a procedure of ultra-violet 
regularization which was implied in Eqs.(|i~6|) and (|i~9p. We start with the NSVZ low 
energy theorem ||23|| , Eq . (|16D , which is here repeated for convenience: 



lim i J dxe^(0\T{^G\x) ^G 2 (0)}|0> = "4<^G 2 > , (A.l) 

where (3{a s ) = — 6a 2 /(27r) + 0(a^) is the Gell-Mann - Low /3-function for YM theory with 
b = (11/3)AT C , and N c is the number of colors. 

The low energy theorem was obtained in [E3l] using the one-loop /3-function with a 



particular attention to a regularization of ultra-violet (UV) divergences, arising in the 
two-point function ( |A.1| ), within Pauli-Villars procedure. In this derivation Dyson type 
of the T-product symbol was implied in Eq.( [A.l[) . It was shown that quadratically diver- 



gent UV contributions cancel out identically in both sides of Eq.( |A.l[) . This implies that 
perturbative contributions should always be subtracted in vacuum condensates such as 
(f3(a s ) / (Aa> s )G 2 ) in Eq. (|A.l|) . Once this rule is accepted, the dependence of any (nonper- 
turbative) condensate (O)np of dimension d on the bare coupling constant go (normalized 
at the cut-off scale Mr) is fixed by the dimensional transmutation formula: 



(O)np = const 



M R exp 



•, 1 \ 1 d 



(A.2) 



and the derivation of the NSVZ theorem proceeds as described in Sect. 3, where the path 
integral definition of correlation functions is used. The latter implies Wick type of the 
T-product symbol. This definition of zero momentum correlation functions flT%D automat- 
ically ensures the same type of renormalization for all such functions, which is fixed by 
a rule of subtracting perturbative UV divergent contributions to the conformal anomaly 
(with e.g. Pauli-Villars regularization). The latter procedure thus defines a nonpertur- 
bative gluon condensate in the ^-vacuum (g 2 G 2 )e, i.e. a nonperturbative part of the 
conformal anomaly calculated for given 9. Its dependence on 1/g 2 is given by Eq.( [A.2| ). 



Zero momentum correlation functions of the operator g 2 G 2 are obtained by the differen- 
tiation of the nonperturbative part of the partition function \og(Z/Z PT ) ( Z PT stands for 
a perturbatively defined partition function which does not depend on 9), where 

Z(0) = Z PT exp{-iVE v (9)} = Z PT exp |— iV<0| - ^G 2 |0) e | , (A.3) 

in respect to \ jg\. Note that the factor d log Z PT /d (g^ 2 ) corresponds to the correlation 
function ( A.l]) in perturbation theory. 



It is a subtraction of perturbative contributions in vacuum condensates (|A.2 ) that we 



here would like to comment upon. Technically, this prescription can be thought of as the 
requirement of absence of regular powers of the coupling constant a s in vacuum conden- 
sates (O)tvp to any finite order in a s . For two-dimensional models, it has been shown 
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that the definition of vacuum condensates via the path integral automatically nulli- 
fies perturbative contributions to the condensates. Moreover, this procedure gives results 
identical to a point-splitting regularization with Dyson type of the T-product. One can 
notice that in four dimensions the separation of genuinely perturbative and nonperturba- 
tive contributions to physical quantities is ambiguous as it depends on a definition of the 
sum of perturbative series. Still, there is nothing wrong with the requirement that the 
"nonperturbative" condensates contain the coupling constant g$ as in Eq.( [A.2|) , while the 



regular powers of g$ are absent. The difference between different regularization schemes is 
reduced in this case to possible finite renormalizations (different numerical values) of the 
vacuum condensates. Such a choice of the nonperturbative gluon condensate in Eq. ( |A.3| ) 
is the only ambiguity for all multi-point correlation functions 

Next we would like to discuss zero momentum correlation functions of the topological 
density operator. With Wick definition of the T-product we obtain from Eq.( 



Dm ij^m^GG^GGm} |0) = = („>, . (A.4) 

(Eq.( |A.4D is written for Minkowsky space, and a stands for the trace of the energy- 
momentum tensor.) A few comments on Eq.( [A.4|) are in order. We note that the two 
definition (through Dyson or Wick T-products) are equivalent for the correlation function 
( |A.1| ). For correlation functions of the topological density Q, this is no longer the case. 
Witten-Veneziano construction || specifically implies |35| Wick T-product in the two- 
point function ( |A.4| ). This can be seen both from the definition of zero momentum 



correlation functions of the operator Q used in ||, and from the fact the a gauge non- 
invariant axial ghost (Veneziano ghost pole) can not appear in Dyson T-product which 
is related to contributions of intermediate gauge-invariant states only. This is why Wick 
type of the T-product was used in Eq.( A.4 ). With this definition the two-point function 



( |A.4| ) does not contain UV divergences which are present in Zpr and drop out after the 
differentiating of log Z in respect to 9. We note that an attempt to calculate the correlation 
function ( |A.4| ) using Dyson T-product (and adding corresponding contact terms) would 
face a problem due to the fact that, in contrast to the trace of the energy-momentum 
tensor, the topological density operator in pure YM theory is not seen to be related 



to any quantity conserved at the classical level. Therefore, the canonical methods with 
Pauli-Villars regularization procedure used in []23| would apparently be not applicable in 
this case. 

We further require that the 9 dependence in Eq.( [A.3|) is described by a single dimen- 
sionless function f{9) such that 

w {a)e = w {{a)o m = {a)o r{6) • (A - 5) 

Any other form of introducing the 9 dependence can be reduced to Eq.( |A.5| ). For example, 
the ansatz (a) s = (a) fi(9) + Ay M / 2 (^) could be transformed to the form ( |A.5| ) by the 
redefinition f\{9) + / 2 (^)Ay M / (cr)o — > f{9). A function f{9) should satisfy the constraints 
f(0) = 1 , f'{0) = 0, which means that its small 9 expansion reads 

/(0) = l-2£ 2 2 + ... , (A.6) 
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where £ 2 is some dimensionless number. Using this in Eq.( |A.4| ), we obtain 



lim i / dxe tqx (0\T ^GG(x) ^GG(0) |O) 0=O = f (- 
<i^o J [on on 



■in 



s G 2 )e=o 



(A.7) 



The assumption made in Eq.([H]) in the text was that Eq.( |A.7p is covariant in 9, i.e. 
remains of the same form ( |A.7| ) not only for 9 = 0, but also for small values 9^0. This 
requirement fixes the function f(9) completely: 



f"{9) = -4ef(0) => f(0) = cos(2£0) 



(A.8) 



which results in the 9 dependence of the vacuum energy and topological density condensate 
identical to the one displayed in Eqs. (^0|) and (^l|). This assumption is self-consistent 



because the same Eqs.(pO|),(pTD can be obtained without fixing the function f{9), but 
instead by postulating the covariant relation (|i~9l) and resumming Taylor expansions in 9 
for the condensates (G 2 )e , (GG)g, as was done in plf . One more relation needed for this 
purpose reads 



dx(0\T{^G 2 ' 



(A.9) 



which is a particular version of the original NSVZ theorem | 23| . Eq.( |A.9|) is valid for any 
small 9. 

Any other choice for the function f(9), different from Eq.( |A.8| ), would presumably not 
be self-consistent, though we did not find a general proof of this statement. One should 
stress that the form ( |A.8| ), implying the reiteration of the parameter £ 2 in multi-point 
correlation functions of the operator Q, is consistent with Veneziano construction H for 
the ghost pole mechanism, where it was found that 



2n-l 



89 



2»- 



-(Q) 



8=0 



>-l)n~l N( 




2n-l 



Nr. 



oo 



A = 0(N° C ) 



(A.10) 



At the same time, Eq.( [A.8|) relates the ghost pole residue in Veneziano scheme with 
conformal anomaly in the theory with 9 ^ 0, in terms of the parameters £ , (cr) . In 
our work [p4| , a particular choice, £ = 4/(36), was advocated using a one-loop connection 
between the conformal and axial anomalies in the theory with an auxiliary heavy fermion. 
Moreover, the covariance of Eq. (|A.7|) in respect to 9 followed automatically within a 
procedure used in pl| , where Eq. (|A.7| ) with £ = 4/(36) was obtained directly from the 
NSVZ relation ( |A.1| ) which is valid for any small 9. 
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